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Spectral Signatures of the Diffusional Anomaly in Water 
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Analysis of power spectrum profiles for various tagged particle quantities in bulk SPC/E water is 
used to demonstrate that variations in mobility associated with the diffusional anomaly are mirrored 
in the exponent of the 1/ f a region. Monitoring of f / f a behaviour is shown to be a simple and direct 
method for linking phenomena on three distinctive length and time scales: the local molecular 
environment, hydrogen bond network reorganisations and the diffusivity. The results indicate that 
experimental studies of supercooled water to probe the density dependence of 1//" spectral features, 
or equivalent stretched exponential behaviour in time-correlation functions, will be of interest. 
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Supercooled water, unlike most liquids, shows an in- 
crease in molecular mobility with compression over cer- 
tain ranges of temperature and pressure 0, 0> 0> 0- 
Along an isotherm, this anomalous diffusional regime 
is bounded by densities correspondin g to the diffusivity 
minimum and maximum [5|, la, 13. la? ISl liil • In this work, 
the diffusional anomaly is related to dynamical varia- 
tions in the special structural features of water: a strong 
preference for local tetrahedral order and the presence 
of a fluctuating three-dimensional network of hydrogen 
bonds. 

We use power spectral analysis to study temporal fluc- 
tuations in various observables. The power spectrum of a 
time-dependent mechanical quantity, A(t), is defined as 



S(f) 



(A(t) (A))e 



2nift 



dt 



(1) 



where (A) is the corresponding average over the system 
trajectory. In the case of water, power spectra associated 
with several quantities show a distinctive 1 / f a type de- 
pendence on the frequency / [111 E1E1E1- Such l/f a 
or flicker noise is a generic feature of systems with multi- 
ple time scales 01 an d, in this context, originates from 
hydrogen bond network rearrangements involving vari- 
able numbers of molecules. This study shows that the 
power spectrum profile associated with various quantities 
sensitive to the local molecular environment, particularly 
in the l/f a regime, carries useful information on spa- 
tiotemporal correlations in the hydrogen bond network, 
including distinct signatures of the diffusional anomaly. 

We have used the extended simple point charge 
(SPC/E) intermolecular potential for water [lfjj since the 
density, diffusional and compressibili ty a nomalies of this 
model are well-characterised HIS El- Molecular dy- 
namics (MD) simulations of bulk SPC/E water in the 
NVT ensemble we re p erformed using the the DL_POLY 
software packa ge Il7l Il8l. Simulation details are sum- 
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Figure 1 shows the power spectra associated with dif- 
ferent tagged particle quantities for bulk SPC/E water at 
230K and 0.9 g cm~ 3 . The tagged molecule potential en- 
ergy is the interaction energy of a given molecule with all 
other molecules in the system. The corresponding power 
spectrum, S u (f), shows a broad peak, centred around 
500 cm -1 , due to the high-frequency, essentially single- 
molecule, librational modes. The S u (f) curve shows two 
regions of l// Q -type behaviour: (i) the high frequency 
60-298 cm" 1 range with a' u = 1.56 ±0.02 and (ii) the low 
frequency 1-40 cm" 1 region with a u — 1.06 ± 0.02. The 
0-0 stretch and 0-0-0 bending modes for bulk SPC/E 
water are known to occur at 200 and 50 cm -1 respec- 
tively [2(J and therefore the high frequency 1/ f a regime 
involves two or three-molecule hydrogen-bond network 
rearrangements. Since increasing delocalisation of vibra- 
tional modes results in lowering of the frequencies, the 
low frequency 1// Q regime must involve displacements of 
four to six molecules. Crossover to Markovian or white 
noise behaviour should occur below 1 cm -1 jlflj which 
is a frequency regime not studied here. Also shown in 
Figure 1 is Sub(f) obtained from fluctuations in the 
librational kinetic energy, defined as the difference be- 
tween the total and the centre-of-mass kinetic energies of 
a rigid molecule. Sub(f) shows no evidence of multiple 
time scale behaviour and the single prominent peak is 
again due to the librational modes. 

The local structural order around a tagged oxygen 
atom is gauged by two order parameters, qo and qu- 
The degree of tetrahedrality of the four nearest hydrogen 
atoms surrounding a given oxygen atom i is measured by: 



= 1 -gE E (cos^ fe + l/3) 2 (2) 
j=i k=j+i 



where tpjk is the angle between the bond vectors ry and 
Yik where j and k label the four nearest hydrogen atoms. 
The qo order parameter is similarly defined using the 
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previously used to characterise orientational order in wa- 
ter [5j. The corresponding power spectra are labelled 
So(f) and Sh(J)- The Sn(f) curve shows a clear 
peak due to librational motion, unlike the So{f) curve. 
Both So(f) and Sh(J) spectra show two 1// Q regimes 
with different exponents but in the case of So{f), the 
higher frequency 1//" regime stretches between 200-1000 
cm . The spectrum of fluctuations in the number of 
nearest oxygen atoms surrounding a tagged oxygen atom, 
Snn(/), is very similar to S (f) @ 

The key features of the power spectral profiles are: (i) a 
broad librational peak (ii) a high-frequency 1/ f a regime 
above 50 cm -1 with exponent a' and (iii) a second, low- 
frequency multiple time scale regime in the range 1- 40 
cm -1 with a different exponent a. A white noise region 
(a = 0), seen at very low frequencies, becomes more pro- 
nounced at higher temperatures and/or higher densities, 
as discussed below. The librational peak is seen in S u (f), 
Su(f) and Sub(f) while the multiple time-scale regions 
are seen in all power spectra except Sub(f). At a given 
state point, the exact values of the exponents in the 1/ f a 
regimes as well as the crossover frequencies depend on the 
specific quantities studied. 

Figure 2 shows S u (f) at different densities for the 
230K isotherm. At this temperature, SPC/E water has a 
diffusional minimum and maximum atp m m=0.9 g cm~ 3 
and pmax—^-^ g cm~ 3 respectively 0,13. The librational 
region of the S u (f) curves changes qualitatively with 
density along this isotherm. It is most clearly demar- 
cated at Pmin=0.9 g cm -3 indicating that decoupling of 
librational modes from the hydrogen-bonded network re- 
arrangements is maximal at this "ice" density. The libra- 
tional peak moves to lower frequencies with increasing 
density and, for p > p ma x, becomes a shoulder. We thus 
provide direct dynamical evidence that the anomalous 
diffusional regime, for which dD/dp > 0, is associated 
with increasing coupling of the librational modes with 
vibrations of the hydrogen-bonded network. This has 
been deduced in earlier work based on static structural 
distributions The density-dependent changes in the 
1// Q regions of the power spectra are quantified by ex- 
amining the associated exponents. Figure 3 shows that 
at 230K, the diffusivity, D, is strongly correlated with the 
low- frequency (1-40 cm -1 ) exponent, a u , and anticorre- 
lated with the high frequency (60-298 cm -1 ) exponent 

To understand the trends in the exponents of the 1/ f a 
regions, it should be noted that power spectra generated 
using a uniform distribution of time scales between fre- 
quency limits Ai and A2, show a 1// regime (a ks 1) for 
Ai < / < A 2 , a white noise regime for / < Ai (a = 0) 
and a Lorentzian tail (a w 2) for / S> A 2 15] . Thus a u 
and a' u will index the shapes of the multiple-time scale 
regime and the Lorentzian tail respectively. A separation 
of two decades in time scales between Ai and A2 is ade- 

nnatp fnr nhsprvi-ncr a 1 / f a rpcrimp with rv plnsp tn 1 At. 



Pmin, the librational modes will be most effectively de- 
coupled from the network vibrations, the range of avail- 
able frequencies contributing to the 1// regime will be 
the narrowest and the Lorentzian tail will be most pro- 
nounced. This is consistent with our observation that a' u 
is maximum and a u is « 1 at p m in- In physical terms, it 
implies that at 230K and close to p m in , vibrational modes 
involving three or less molecules are significantly decou- 
pled from more delocalised network re-organisations. We 
focus on the correlation with the diffusivity of the ex- 
ponent a u associated with the low-frequency, multiple 
scale region of the power spectrum. As librational modes 
become progressively more coupled to the network vibra- 
tions, high frequency components are introduced result- 
ing in an increase in a u with density in the anomalous 
regime between p m in and p ma x ■ The decrease in a u with 
increasing density for p > p max is due to increasing im- 
portance of steric effects which force the system to behave 
more like a simple liquid with a = (l2T |. 

The behaviour of the librational peak in Sub(f) and 
Sn(f) at 230K is similar to that seen for S u (f) and the 
density dependence of the exponents in the l// Q regions 
is shown in Figure 3. The olh and oto values are close 
to 1 in both cases. The density-dependent variations in 
these quantities, an and ao, are smaller in magnitude 
than those of a u though they show similar trends. 

Increasing temperature or density results in coupling 
of librational modes to the network vibrations and at 
300K, the librational band profile is virtually invariant 
with density. At 260K, two 1// Q regimes with different 
slopes can be seen with a low-frequency regime between 
1-40 cm -1 and a high-frequency Lorentzian tail between 
60-200 cm -1 . At 300K, only a single 1// Q regime can be 
seen between 4 and 200 cm -1 with a crossover to white 
noise behaviour between 0.1 and 0.5 cm -1 . Figure 4 il- 
lustrates very clearly the correlation between the density- 
dependent variations in the diffusivity and the exponent 
a u . The attenuation of the diffusional anomaly with tem- 
perature and its disappearance at 300K are reflected in 
the behaviour of a u . 

We have also examined the So(f) and <5f/(/) spectra 
at 260K and 300K. Compared to S u (f), the crossover 
to white noise behaviour occurs at higher frequencies for 
these quantities; for example, between 5 to 10 cm -1 , for 
So{f) and Sn(f) at 300K. This indicates that local 
tetrahedral order decorrelates faster than fluctuations in 
the tagged molecule potential energy. 

This study establishes that the variations in mobility 
with pressure in the region of the diffusional anomaly are 
unambiguously mirrored in the exponent of the 1//" re- 
gion thus allowing one to monitor the dynamical effects of 
coupling the localised librational modes to the hydrogen- 
bonded network. The range of frequency over which 
1// Q behaviour is observed, including the frequency of 
crossover to white noise, provides a quantitative assess- 
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tions in the network vanish. Moreover, l/f a behaviour is 
not equally pronounced for all observables, due to differ- 
ing sensitivities to processes operating on different time 
scales. 

The significance of the present results derive from the 
clear connections that can be drawn between between 
three distinctive length and time scale features of bulk 
water: the local molecular environment, the hydrogen 
bond network reorganisations and the diffusivity which 
is a long-time averaged transport property. This ap- 
proach provides a possible route for understanding the 
connection between structural order and the static and 
dynamic anomalies of water 0, l2l) and can be extended 
to other hydrogen-bonded systems or networked liquids. 
Our work also suggests that it would be of interest to ex- 
amine the connection between 1/ f a behaviour and spa- 
tiotemporal heterogeneity in supercooled liquids since the 
distribution of length and time scales in the system is 
expected to widen as one approaches the kinetic glass 
transition. In the context of water, such a kinetic glass 
transition is predicted for SPC/E water around 190K at 
1 g cm " 1 jlSlifl; the experimental value is around 225K 
[2J, |23j . Near this dynamical transition, water behaves in 
a very similar manner to the binary Lennard-Jones glass- 
formers with no apparent effect of the hydrogen bond 
network. Our I// Q analysis shows that by 230K, two 
and three-molecule vibrational modes are shifted into the 
Lorentzian tail, indicating that these vibrational frequen- 
cies are progressively decoupled from overall network re- 
organisations. As the temperature is reduced further, 
one would expect the multiple time-scale regime to shift 
to lower frequencies and involve larger clusters; the dy- 
namics of these larger clusters may be insensitive to di- 
rectional bonding and local tetrahedral order. 

Our results have interesting implications for experi- 
mental work on water and aqueous solutions, specially 
in the context of hydration of biomolecules [24], |2jj • We 
show that different observables will have different degrees 
of sensitivity to the underlying multiple time-scale dy- 
namics of the hydrogen-bonded network and therefore 
the choice of spectroscopic technique is likely to be sig- 
nificant. For example, the behaviour of high-frequency 
librational modes as a function of density can be stud- 
ied by ultrafast spectroscopy . Experimental evidence 
for 1//" behaviour, or the equivalent stretched exponen- 
tial behaviour of time-correlation functions, exists from 
inelastic neutron scattering, Raman and, most recently, 
optical Kerr effect measurements [H Hi El 0. The 
pressure or density dependence of the exponents has not 
been examined so far and our work indicates that these 
will be of considerable interest . 
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Figure Captions 



1. Power spectra associated with temporal fluctu- 
ations in different tagged particle quantities at 
230K and 0.9 g cm" 3 . The different curves rep- 
resent S u (/) (thick solid line), Sm(f) (solid line), 
Snn(J) (dot-dashed line), So(f) (dotted line) and 
Siiif) (dashed line). Arrows show the 0-0 stretch 
(200 cm" 1 ), 0-0-0 bend (50 cm" 1 ) and the libra- 
tional region (400-1000 cm -1 ) 2Q|. Systems at a 
density of 0.75g cm -3 are likely to be inhomoge- 
neous Q- 



2. Power spectra, S u (f), at different densities (in g 
cm" 3 ) along the 230K isotherm. 
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3. Behaviour of the exponent associated with the 1// Q 
regimes as a function of density at 230K. Part 

(a) shows the exponents associated with the high- 
frequency 1// Q regime with a' u , a' H and a' Q eval- 
uated over the frequency ranges 60-298 cm -1 , 90- 
298 cm -1 and 200-1000 cm" 1 respectively. Part 

(b) shows the exponents associated with the low- 
frequency 1/ f a regime with a u , an and ao evalu- 
ated over the the 1-40 cm" 1 range. Also shown is 
the diffusivity, D, in units of 2 x 10 -6 cm 2 /s. 



4. Dependence on density, p, of (a) the diffusivity, D, 
and (b) the exponent a u along isotherms at 230K, 
260K and 300K. The exponent a u is evaluated us- 
ing the low- frequency 1/ f a regime from 1-40 cm -1 
at 230K to 280K and over the frequency range 4- 
200 cm" 1 at 300K. 
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